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Abstract 

As part of a general study of the time-dependent local density approximation 
(TDLDA), we here report calculations of optical activity of chiral molecules. 
The theory automatically satisfies sum rules and the Kramers-Kronig rela- 
tion between circular dichroism and optical rotatory power. We find that 
the theory describes the measured circular dichroism of the lowest states in 
methyloxirane with an accuracy of about a factor of two. In the chiral fullerene 
C76 the TDLDA provides a consistent description of the optical absorption 
spectrum, the circular dichroism spectrum, and the optical rotatory power, 
except for an overall shift of the theoretical spectrum. 
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I. INTRODUCTION 



The time- dependent local density approximation (TDLDA) as well as the time-dependent 
Hartree-Fock theory has been applied to the optical absorption of atomic and molecular 
systems with considerable success [1-14]. Here we want to see how well the TDLDA method 
does on a more subtle aspect of the optical response, the optical activity of chiral molecules. 
Calculation of circular dichroism and especially optical rotatory power is more challenging 
because the operators that must be evaluated are more sensitive to the approximations on 
the wave function than the electric dipole in the usual form ef. Nevertheless, we anticipate 
that TDLDA may be a useful theory because the operators are still of the form of single- 
electron operators. The TDLDA is derived by optimizing a wave function constructed from 
(time-dependent) single-particle wave functions, so its domain of validity is the one-particle 
observables. 

In our implementation of TDLDA ||, we represent the electron wave function on a uni- 
form spatial grid. The real-time evolution of the wave function is directly calculated and 
the response functions are calculated by the time-frequency Fourier transformation. The 
method respects sum rules and the Kramers-Kronig relation between the circular dichroism 
and optical rotatory power. Since the grid representation is bias-free with respect to elec- 
tromagnetic gauge, it is not subject to the gauge difficulties encountered when the space of 
the wave function is constructed from an atomic orbital representation. 

Optical activity has been a challenging problem for computational chemistry, but there 
has been considerable progress in recent years. For example, Carnell et al. (T^j present a good 
description of the circular dichroism of excited states of R-methyloxirane using a standard 
Gaussian representation of the wave function. The optical rotatory power is a much more 
difficult observable, since the whole spectrum contributes. Only very recently have ab initio 



calculations been reported for this property [16 , 17] . 

After presenting our calculational method, we report our exploratory study on optical 
activities of two chiral molecules: R-methyloxirane, a simple 10-atom molecule with known 
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chiroptical properties up to the first few excited states [|T^,|T^,|T^j , and C76, a fullerene with 
very large optical rotatory power and significant circular dichroism in the visible and UV 
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II. FORMALISM 



A. Some definitions 



Polarization of chiral molecule in applied electromagnetic field is expressed using two 



coefficients a and (3 as 21 



P=aE-?-. (1) 
Here a is the usual polarizability and is given microscopically as 

where $ n and _E n are the eigenvector and eigenvalue of the ra-th eigenstate of the many- 
body Hamiltonian H, H$ n = E n $> n , and E n Q = E n — Eq. The 5 is an infinitesimal positive 
quantity. Employing the oscillator strength 

fn= 2 -^ 1 -(%\E^n) 2 , (3) 

we define the optical absorption strength whose integral is normalized to the active electron 
number, 

S(E) = ^5(E-E n )f n . (4) 

n 

It is related to the imaginary part of the polarizability, 



= ^Imaf^) 

n, e 2 7r 



The basic quantity which characterizes the chiroptical transition is the rotational strength 
defined by [^2 



R n = -pL { $ \J2n\® n ) ■ ($ B |£fx V|$ ). (6) 



2mc 

We define the complex rotational strength function, 

K{E) = ¥ Uo-S-i* " J ^ (7) 

The beta function in eq.(l) is related to TZ(E) by /3(E) = Mj1Z(E). We will also use the 
rotational strength function R(E) defined by 

R(E) = E n )R n = ^iS. (8) 

7T 

n 

As is seen below, the optical rotatory power is proportional to the real part of (3(E), and the 
circular dichroism to R(E). They are related to each other by the Moscowitz's generalized 
Kramers-Kronig relation [p3[|. 

The difference of complex indices of refraction for left and right circularly polarized light 
is proportional to 1Z in dilute media; the relation is 

n L -n R = ^TZ(E), (9) 

where Ni is the number of molecules per unit volume. For comparison with experiment, the 
common measure of circular dichroism is the decadic extinction coefficient, given by 

Ae= A cm Clog e 10 ImK -^ ) ' (10) 

where C is the concentration of molecules in moles/liter and the subscript on the wavelength 
A is a reminder to express it in centimeters. The optical rotatory power is conventionally 
reported as 

[a] = 180°-^-Re(n L - n R ), (11) 
where C gm is the concentration of molecules in gm/cm 3 . 
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B. Real-time TDLDA 



We first rewrite the above strength functions as time integrations. We employ the 
time dependent wave function ty(t) = exp[— iHt/h]^(0) with the initial wave function at 
t — given by ^(0) = exp[ik J2i z i\®o, where A; is a small wave number. In the linear 
response, the time-dependent polarizability is proportional to the dipole matrix element, 
z(t) = (&(t)\ J2i z i\^(t))- The frequency dependent polarizability in z direction is then 
obtained as the time-frequency Fourier transformation of z(t), 

a z {E) = ^ r ^ E+i5 ^ h z(t). (12) 
k Jo n 

The polarizability a(E) is given by the orientational average, a = (a x + a y + a z )/3. 

Similarly, we denote the angular momentum expectation value as L z (t) = ($(t)\ — i(rx 
V) z \^(t)). To linear order in k, we may express it as 

L z (t) = -2^ COS (^f) ($ |E^I $ n)($n|E(^ >< ^i),|*0>- (13) 

The complex rotatory strength function 71(E) is expressed as, 

^Me^W^-E^ETJs) <*oII>I*»><*„IX>->< V).l*o>. (14) 
71(E) of eq.(7) is the sum over three directions, TZ = 7Z X + 7Z y + 7Z Z . 

In the time-dependent local density approximation, the time-dependent wave function 
ty(t) is approximated by a single Slater determinant. We prepare the initial single electron 
orbitals as ^(O) = exp[ikz](pi where 4>i is the static Kohn-Sham orbitals in the ground state. 
The ipi(t) follows the time- dependent Kohn-Sham equation, 

ih^i(t) = {-^ V2 + E V ^n(r~ Ra) + e 2 j dft^L + » xc (p(f, t)) } tk(t), (15) 

where Vi on is the electron-ion potential and /i xc is the exchange-correlation potential. The 
time-dependent density is given by p(f,t) = \^(t, t)\ 2 . The time-dependent dipole mo- 
ment may be evaluated as z(t) = X^V'iMV'i) an d the similar expression for L z (t). The 
strength functions are then evaluated with eqs.(12) and (14). 
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C. Sum rules 



According to the TRK sum rule, the integral of S(E) is equal to the number of active 
electrons N. This sum rule is respected by the TDLDA. It also appears in the short time 
behavior of z(t) as 

hk 

z (t) = N — t (t small) (16) 
m 

For the rotational strength, we define energy-weighted sums as 

RW = r dEE n R(E). (17) 
J o 

It is known that R^ for n < 4 vanishes identically in the exact dynamics [22| , 24 -|26|. The 



vanishing of R^ in the time-dependent Hartree-Fock theory was first noticed in [p7 ]. The 
short time behavior of L(t) = L x {t) + L y (t) + L z (t) is related to R^ n ' as 

Here we note that L(t) is an even function of t as seen in eq.(13). Since R^ = for n < 4, 
we see that L(t) behaves as t 6 for small time. Li(t)(i = x,y,z) behave as t 2 at small t and 
the cancelation up to t A order occurs after summing up three directions. In the TDLDA 
dynamics, we confirmed that at least t° and t 2 coefficients of L(t), namely R^ and R^ 2 \ 
vanish identically. 

D. Numerical detail 

The TDLDA uses the same Kohn-Sham Hamiltonian as is used in ordinary static LDA 
calculations. As is common in condensed matter theory, we use pseudopotentials that in- 
clude the effects of i^-shell electrons rather than include these electrons explicitly. The 
pseudopotentials are calculated by the prescriptions of ref. |28[] and |29[. We employ the 



simple exchange-correlation term proposed in ref. pD|rjT| ]. There are improved terms now in 
use [|3|,[3^], but it was not deemed important for our application. 
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There are many numerical methods to solve the equations of TDLDA. Ours uses a 
Cartesian coordinate space mesh to represent the electron wave functions, and the time 
evolution is calculated directly |||9],|34] . There are only four important numerical parameters 
in this approach: the mesh spacing, the number of mesh points, the time step in the time 
integration algorithm, and the total length of time that is integrated. We have previously 
found that the carbon molecules can be well treated using a mesh spacing of 0.3 A PjMf. 
We find 0.25 A is necessary for methyloxirane to represent accurately the orbitals around 
oxygen. We take the spatial volume to be a sphere of radius 8 A both for methyloxirane 
and C76 presented below. The total number of mesh points, defining the size of the vector 
space, is about 4vr J R 3 /(3Ax 3 ) ~ 80,000(140,000) for mesh size of 0.3 A(0.25 A). 

The algorithm for the time evolution is quite stable as long as the time step At satisfies 
At < h/\H\, where \H\ is the maximum eigenvalue of the Hamiltonian. This is mainly 
dependent on the mesh size. For Ax = 0.3 A, we find that At = 0.002/i/eV is adequate. 
We integrate the equation of motion for a length of time T = 60h/eV for C76 (50/i/eV for 
methyloxirane). Then individual states can be resolved if their energy separation satisfies 
AE > 2nh/T ~ O.leV. 

Our numerical implementation, grid representation of the wave function and the time- 
frequency Fourier transformation for the response calculation, has several advantages over 
the usual approach using basis functions centered at the ions. They include: 

(1) The full spectrum of wide energy region may be calculated at once and it respects sum 
rules. The non-locality of the pseudopotential may cause violation of the sum rule, but the 
effect is small in the present systems. 

(2) Since the circular dichroism R(E) and the optical rotatory power, real part of (3(E), are 
calculated as Fourier transformation of single function L(t), the Kramers-Kronig relation is 
automatically satisfied. 

(3) The gauge independence of the results is satisfied to high accuracy. Employing the 
commutation relation [H,J2i^i] — ~~ IrX^Vi, there is alternative expressions for optical 
transitions with gradient operator instead of coordinate. For the rotational strength, for 
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example, 

R n = - o ( $ 0lEV^n) ■ ($n|E^ X Vi|$0>, (19) 

zm ch, n Q i i 

The strength function with this expression may be calculable with initial wave function 
ipi(0) = exp[idV z ]<f>i with small displacement parameter d. Since the grid representation of 
wave function does not have any preference on the gauge condition, our method gives almost 
identical results for the coordinate and gradient expressions of dipole matrix elements. 



III. R-METHYLOXIRANE 



The geometry of R-methyloxirane is shown in Fig. 1. We use the same nuclear distances 



as in ref. [f[q| . We show in Fig. 2 the results of the static calculation for the orbital energies. 
We find a LUMO 6.0 eV above the HOMO, and a triplet of unoccupied states 0.5 eV higher. 
In our calculation the lowest unoccupied orbitals have an diffuse, Rydberg-like character, s- 
wave for the lower and p-wave for the upper, as in previous calculations [|I5| , |35|| . The HOMO 
is localized in the vicinity of the oxygen atom, and the measured absorption strength seen 
at 7.1 and 7.7 eV is attributed to the excitation of a HOMO electron to the diffuse states. In 
the TDLDA, the excitation energy comes out close to the orbital difference energies, except 
for strongly collective states. Indeed we find in the TDLDA calculation that the excitations 
are within 0.1 eV of the HOMO-LUMO energy and the energy difference for the next state 
above the LUMO. This is one eV less than the experimental values. It is known that the 
LDA energy functional that we use is subject to overbind excitations close to the ionization 



threshold. There are improved energy functions that rectify this problem ||36|| , but for this 
work we judged the error not important. 

The next property we examine is related to the electric dipole matrix element, namely 
the oscillator strength / associated with the transition. The optical absorption strength is 
shown in Fig. 3. The total strength up to 100 eV excitation is f=22.4, which is 93% of the 
sum rule for the 24 active electrons. Notice the lowest two peaks, centered at 6.0 and 6.5 
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eV. These are the states we are interested in. Their oscillator strengths are given in Table 
I. We see that the states are both weak, less than a tenth of a unit. The effect of the time- 
dependent treatment is to lower the strengths by 30-50%. This is the well-known screening 
effect associated with virtual transitions of more deeply bound orbitals. We find that the 
computed transition strengths are within a factor of two of the measured ones. Typically, 
the TDLDA does somewhat better than this, but most of the experience has been with 
transitions carrying at least a tenth of a unit of oscillator strength. The original theoretical 



calculation gave very poor transition strengths fl8fl , off by more than an order of magnitude. 



Unfortunately, the more recent study |L5j did not include theoretical transition strengths. 



We numerically confirmed that our method gives almost identical results with coordinate 
and gradient expressions of dipole matrix elements, as we noticed in the previous section. 
However, exceptionally, the oscillator strength of the very weak features discussed above 
suffer substantial dependence on the expression. With gradient formula for the transition 
matrix elements, strengths of both first and second transitions are larger by about factor 
two than the coordinate expression. Since the gradient formula emphasizes high-momentum 
components more heavily, we think the results with coordinate matrix elements may be more 
reliable for low transitions, and we quote them in Table I. 

We now turn to the chiroptical response. Fig. 4 shows the short-time behavior of L x (t) 
and the sum of the three Cartesian components L(t) = J2{Li(t). An initial perturbation 
of k = 0.001 A -1 is employed. To within numerical precision, L x (t) (solid line) grows with 
time as t 2 , as discussed below eq.(18). The same is true for the other two components, L y 
and L z . This shows that the numerical algorithm respects the first sum rule. The combined 
response, L(t) (dashed lines) shows an extreme cancelation at short times, as required by 
the additional sum rules. However our numerical accuracy does not allow us to determine 
the order of the power behavior. The evolution of L(t) for larger times is shown in Fig. 5. 
Physically, the TDLDA breaks down at long times because of coupling to other degrees of 
freedom. A typical width associates with such couplings is of the order of a tenth of an eV, 
implying that the responses damp out on a time scale of T rs lOh/eV. We note that the 
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TDLDA algorithm itself is very stable, and allows us to integrate to much larger times and 
obtain very sharp theoretical spectra. 

We next show the Fourier transform of the chiroptical response. The circular dichro- 
ism spectrum calculated with eq. (14) is shown in Fig. 6. Here we have integrated L to 
T = 50h/eV, including a filter function in the integration to smooth the peaks. One can see 
that the s- and the p-transitions are clearly resolved, although the three p-transitions are not 
resolved from each other (as is the experimental case). The s-transition has a negative circu- 
lar dichroism and the p-transition a positive one. Integrating over the peaks, the strengths 
of the two peaks are -0.0014 and +0.0014 A 3 -eV, respectively. The strengths are commonly 
quoted in cgs units; the conversion factor is 1 eV-A 3 = 1.609 xlO" 36 erg-cm . The values in 
cgs units are given in Table I, compared to experiment and previous calculations. We find 
the signs are correctly given, but the values are somewhat too high, by a factor of 2 or 3. 
The calculation of ref. Jl5| gave a result within the experimental range for the p-multiplet 



but too small (by a factor of 2) for the s-transition. Thus we find that the TDLDA has a 
somewhat poorer accuracy in this case. 

Next we consider the optical rotatory power. It could be calculated as the real part of the 
Fourier transformation eq.(14). In practice, however, we found the calculation employing 
Kramers-Kronig relation to the rotational strength function, 

Re(3(E) = 2 fcJ o °°dE'^^- 2J (20) 

gives more accurate result especially at the energy below the lowest transition. The mea- 



surement is available at sodium D-line, 2.1 eV, [a]o = +14.65° which gives (3 = +0.0017 
A 4 . The calculated value at low energy is very sensitive to the number of states included in 
the sum. Fig. 6 shows the calculated value as a function of a cutoff energy, upper bound in 
the integration in eq.(20). The value taking only the contribution of the lowest transition 
is -0.06. Including more states produces values that fluctuate in sign and magnitude within 
that range. Including all states below 100 eV gives a cancelation by a factor of 60 to yield 
a value (3 = —0.001 A 4 . This has the opposite sign but the same order of magnitude as the 
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measured (3. Clearly, to get high relative accuracy with such a strong cancelation is more 
demanding than our TDLDA can provide. 



IV. C 



71 i 



Remarkably, it is possible to separate the chiral partners of the double-helical fullerene 
C76 using stereospecific chemistry fl20 |. The molecule shows a huge optical rotatory power, 



[ot\r> = —4000°, and a complex circular dichroism spectrum between 2 and 4 eV excitation 
|20| . There has been reported only semi-empirical quantum chemistry calculation for the 
optical activity of this molecule 13] 



We first remark on the geometry of the molecule, which has a chiral D2 symmetry p8 
The accepted geometry is depicted in ref. p0| ; it may be understood as follows. We start 
with Cgo, in which all carbons on pentagons. Group the pentagons into triangles and divide 
the fullerene in half keeping two adjacent triangles of pentagons intact in each half. The 
"peel" of six pentagons already has a chiral geometry dependent on the relative orientation 
of its two triangles of pentagons. The C76 is constructed by adding 16 carbon atoms between 
the split halves of the C§$. The added carbon atoms lie entirely on hexagons which form 
a complete band around the fullerene. The inserted band has the geometry of an (8,2) 
chiral buckytube. The result is then the chiral C76. Our calculations are performed on a 
right-handed C76, in the sense that the band of hexagons corresponds to a right-handed 
buckytube. This is the same convention as used in ref. ||38|| , their Fig. 3d. 

The C76 has 152 occupied spatial orbitals. We show the orbitals near the Fermi level in 
Fig. 8. The HOMO-LUMO gap is only 0.9 eV, and there are many transition in the optical 
region. In Fig. 9 we show the optical absorption strength function for the range 0-50 eV. 
Smoothing is made with the width of 0.2 eV in the Fourier transformation. A concentration 
of strength is apparent at 6 eV excitation; there is a similar peak in graphite and Ceo which is 
associated with rr — it* transitions. The strong, broad peak centered near 20 eV is associated 
with cr — c* transitions and is also present in C 60 The feature at 13 eV is not present in 
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Ceo , however. In the next figure, Fig. 10, we show a magnified view of the absorption at low 
energy. We also compare the TDLDA strength with the single-electron strength, smoothed 
also by convolution with a Breit-Wigner function of width T = 0.2 eV. The TDLDA has a 
strong influence on the strength distribution, decreasing the total strength in the low energy 
region and concentrating in the 6 eV peak. The experimental absorption strength |5^] (with 
arbitrary normalization) is shown as the dashed line. It agrees with the TDLDA quite well. 

We next examine the circular dichroism spectrum. Fig. 11 shows the rotatory strength 
function between and 50 eV. Like the case of methyloxirane, it is irregular without any large 
scale structures. Its integral is zero to an accuracy of 0.001 eV-A 3 . The low energy region is 



shown in Fig. 12. Here one sees qualitative similarities between theory and experiment [20 



The theoretical sharp negative peak at 1.8 eV corresponds to an experimental peak at 2.2 
eV. Shifting the higher spectra by that amount (0.6 eV), one sees a correspondence between 
the next positive and negative excursions. We note that a similar shift in the excitation 
energy was also seen in the optical absorption of Cqq between the measurement and the 
TDLDA calculation g. 

Our theoretical optical rotatory power is plotted in Fig. 13. The situation here is dif- 
ferent from the methyloxirane, in that rotatory power is large in a region where there are 



allowed transitions. The measured optical rotatory power, [oc\d = —4000° at 2.1 eV |20 
corresponding to f3 = —7.3 A 4 , is shown as the star. It does not agree with theory, but we 
should remember that the spectrum needs to be shifted by 0.6 eV to reproduce the circular 
dichroism. Adjusting the theoretical spectrum upward by that amount, we find that it is 
consistent in sign and order of magnitude with the measurement. Since the optical rotatory 
power in the region of allowed transitions changes rapidly as excitation energy, a measure- 
ment of the energy dependence would be very desirable, and would allow a more rigorous 
comparison with the theory. 
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V. CONCLUDING REMARKS 



We have presented an application of the time-dependent local density approximation 
to the optical activities of chiral molecules. Our method is based on the uniform grid 
representation of the wave function, real-time solution of the time-dependent Kohn-Sham 
equation, and the time-frequency Fourier transformation to obtain the response functions. In 
this way we can calculate the optical absorption, circular dichroism, and the optical rotatory 
power over a wide energy region, respecting sum rules and Kramers- Kronig relation. 

We applied our method to two molecules, methyloxirane and C 76 . For lowest two tran- 
sitions of methyloxirane, the TDLDA reproduces the absorption and rotational strength in 
the accuracy within factor two. The qualitative feature of the circular dichroism spectrum of 
C 76 is also reproduced rather well. However, the optical rotatory power is found to be a very 
sensitive function with strong cancelation. Even though we obtained rotational strength of 
full spectral region, it is still difficult to make quantitative prediction of optical rotatory 
power at low energies in our present approach. 
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TABLES 



TABLE I. Transitions in R-methyloxirane 



Level 


TDLDA LDA(Free) 


Other theory 
ref. @ ref. @ 


Experiment 
ref. | IE] ref. [18] 


1 E (eV) 
/ 

R (xlO 40 cgs) 


6.0 5.97 
0.012 0.021 
-23.0 -10.4 


6.25(7.12) 6.40 
0.0004* 
-6.43 -2.66*'** 


7.07 7.12 
0.025 

-12.56 -11.8** 


2-4 E (eV) 

Si? (xlO 40 cgs) 


6.5 6.55 
0.044 0.069 
23. 29.9 


6.95 7.3 

0.0012* 
7.93 2.24*'** 


7.70 7.75 
0.062 

6.98 10.8** 



* Calculation with coordinate expression of dipole moment. 
** Negative of value for S-methyloxirane. 
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FIGURE CAPTIONS 



Fig. 1. View of R-methyloxirane with hydrogen on the chiral carbon in the back (and 
not seen). The chirality is R because the three other bonds are arrange clockwise in the 
sequence O, CH 3 , CH 2 . 

Fig. 2 Static LDA orbitals in methyloxirane. 

Fig. 3 Optical absorption strength of methyloxirane in the energy region 0-50 eV. 
Fig. 4 Short-time chiroptical response of R-methyloxirane. Solid line is L x (t), dashed 
line is 

Fig. 5 Chiroptical response of R-methyloxirane L(t) for longer times. 
Fig. 6 R in R-methyloxirane in the interval 0-20 eV. 

Fig. 7 Optical rotatory power Re/3 at E — 2.1 eV as a function of cutoff energy E max in 
the integration of eq. (20). 

Fig. 8 Static LDA orbitals in C76 near the Fermi level. 

Fig. 9 Optical absorption spectrum of C76 in the range 0-50 eV. 

Fig. 10 Optical absorption spectrum of C76 in the range 0-8 eV. Dotted line is the 
single-electron strength, solid line the TDLDA, and dashed line experiment |38| . 
Fig. 11 Circular dichroism spectrum of C76. 

Fig. 12 Circular dichroism spectrum of C76 comparing theory (solid line) and experiment 
(dashed line). The experimental data is from ref. [2(J and is with arbitary normalization. 

Fig. 13 Optical rotatory power of C76, given as Re/3 in unit of A 4 . The cross is the 
measured value from [a]z>. 
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FIGURES 




FIG. 1. View of R-methyloxirane with hydrogen on the chiral carbon in the back (and not 
seen). The chirality is R because the three other bonds are arrange clockwise in the sequence O, 
CH 2 0,CH 3 . 
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FIG. 2. Static LDA orbitals in methyloxirane. 
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FIG. 3. Optical absorption strength of methyloxirane in the energy region 0-50 eV. 
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FIG. 4. Short-time chiroptical response of R-methyloxirane. Solid line is L x (t), dashed line is 
J2iLi(t), in units of A. 
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FIG. 5. Chiroptical response of R-methyloxirane L(t) for longer times. 
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FIG. 6. R in R-methyloxirane in the interval 0-20 eV. 
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FIG. 7. Optical rotatory power Re/3 at E = 2.1 eV as a function of cutoff energy E max in the 
integration of eq. (20). 
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FIG. 8. Static LDA orbitals in C76 near the Fermi level. 
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FIG. 9. Optical absorption spectrum of C76 in the range 0-50 eV. 
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FIG. 10. Optical absorption spectrum of C76 in the range 0-8 eV. Dotted line is the sin- 



de-electron strength, solid line the TDLDA, and dashed line experiment [38]. 
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FIG. 12. Circular dichroism spectrum of C76 comparing theory (solid line) and experiment 



(dashed line). The experimental data is from ref. [20] and is with arbitary normalization 
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FIG. 13. Fig. 13 Optical rotatory power of C76, given as Re/3 in unit of A 4 . The cross is the 
measured value from [a]o- 
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